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ABSTRACT 

*-£h A survey of the Galactic plane in the region —60° ^ / ^ 30°, \b\ ^ 0.25° was carried out 

using the seven-beam Parkes Methanol Multibeam (MMB) receiver, which operates 
Q at a frequency of 6.5 GHz. Three pulsars were discovered, and 16 previously known 

%-[ pulsars detected. In this paper we present two previously-unpublished discoveries, 

c/3 both with extremely high dispersion measures, one of which is very close, in angular 

distance, to the Galactic centre. The survey data also contain the first known detection, 
at radio frequencies, of the radio magnetar PSR J1550— 5418. Our survey observation 
was made 46 days prior to that previously published and places constraints on the 
J> beginning of pulsed radio emission from the source. 

£C) The detection of only three previously undiscovered pulsars argues that there 

t" — ■ are few pulsars in the direction of the inner Galaxy whose flux density spectrum 

is governed by a flat power law. However, these pulsars would be likely to remain 
undetected at lower frequencies due to the large amount of scatter broadening which 
affects pulsars with high values of dispersion measure. Surveys with future telescopes 
C^) at high observing frequencies will, therefore, play an important role in the discovery of 

<^> pulsars at the Galactic centre. By simulating pulsar surveys of the Galaxy with Phase 

^ — I 1 SKA at frequencies of 1.4 GHz and 10 GHz, we find that high-frequency observations 

are the only way to discover and observe the Galactic-centre pulsar population. 



X 



Key words: pulsars: general - stars: neutron - pulsars: individual: PSR J1834-0812 
- pulsars: individual: PSR J1746-2850 



1 INTRODUCTION 

Pulsars are known to have steep spectra in the radio regime. 
At frequencies above about 100 MHz, the relationship be- 
tween flux density, S v , and observing frequency, v, can be 
approximated by the power law S v oc v a with a spectral 



index a. From a large sample of pulsars, |Maron et al. (2000) 
derived a mean value for a of —1.8. For this reason, pul- 
sar surveys have frequently been conducted at frequencies 
<1500 MHz - for example the Parkes southern pulsar survey 
( Manchester et al.|[l9^5| at 436 MHz and the G reen Bank 
pulsar survey at 350 MHz ( jHessels et al.| [2008), while the 
higher frequency of 1400 MHz was used in the Swinburne 
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intermediate-latitude pulsar survey (Edwards et al. 20011, 



Manchester 



the Parkes multibeam pulsar survey (PMPS 
et al.|200~ I and the PALFA survey at Arecibo ( Cordes et al 
20061. In addition, the beamwidth of a telescope scales as 



v ' and so the amount of time required to survey a given 
area of sky with a fixed time per pointing is reduced. 

Charged particles in the interstellar medium cause a 
signal of frequency / to be delayed, relative to a signal at 
infinite frequency, by a time 



t 



DM 



27rm e c f 2 ' 

where we define the dispersion measure, DM, to be 

rd 



DM : 



n c dl 



(1) 



(2) 



2 S. D. Bates et al. 



for an electron number density n c and path length from the 
pulsar to Earth, d. At low frequencies, therefore, the dis- 
persive delay is large, making it easier to discern genuine 
astrophysical signals at low DMs from terrestrial radio fre- 
quency interference (RFI). 

However, at low frequencies, there are two factors which 
can have a strong limiting effect on the sensitivity of a sur- 
vey. The first is the v~ 2S dependence of the sky background 
temperature, T s k y , (e.g. Lawson et al. 1987), which arises 



Table 1. Observational parameters for the MMB pulsar survey. 



due mainly to synchrotron radiation from free electrons in 
the Galactic magnetic field. This can dominate the system 
temperature, T sys , at low frequencies, reducing sensitivity 
to all pulsars - especially those at low Galactic latitudes. 

The second of these factors is interstellar scattering, 
which causes a signal to be broadened with a typically ex- 
ponential tail of timescale r s , which scales with frequency 



approximately as r s oc v (Lohmer et al 



2001 



Bhat et al. 



2004). Unlike dispersion, this scatter broadening cannot be 



easily removed in the signal processing. Therefore, at higher 
frequencies a survey of the Galactic plane will be less af- 
fected by smearing of pulses from distant sources - albeit at 
the cost of reduced intrinsic flux density, but mitigated by 
reduced sky background temperature and larger available 
bandwidth. 

This effect is important, and provided two motivations 
for undertaking this high-frequency survey of the Galac- 
tic plane. Firstly, the distribution of pulsars in the inner 
Galaxy is not well known. Following earlier work by | John^] 

( 2006 1 showed that there may 



ston 



(1994), Lorimer et al. 



be a dearth of pulsars in this region, but that this result is 
highly dependent on the poorly-known free-electron distri- 
bution. Finding pulsars in this region, through surveys at 
higher frequencies (e.g. |Johnston et al.|2006[ ), can constrain 
both the distribution of pulsars and dispersive material. 

Secondly, it has long been known that many exciting 
tests of General Relativity, as well as a large amount of evo- 
lutionary information, would be facilitated by the discovery 
of a binary system containing both a pulsar and a black hole 
( Kramer et al.|2004 l. By looking towards the Galactic cen- 
tre, one might hope to observe a pulsar in a long orbit around 
the supermassive black hole, or compact stellar-mass bina- 



ries produced in the high stellar density conditions (Pfahl & 

[Loeb|[2004| . 

In this paper we present the findings of the first wide- 
area pulsar survey at a frequency above 5 GHz. In Section 
2 we discuss the parameters of the survey and the resulting 
sensitivity limit of our observations, and briefly mention the 
processing of the data. In Section 3 we explain the problems 
that arise from searching for pulsars at this high frequency 
and the techniques used to minimise these difficulties. In 
Section 4 we outline the parameters of the newly discov- 
ered pulsars, and then look at redetections of known pulsars 
in the survey area and use these to measure the variation 
of scattering timescale with dispersion measure (hereafter 
DM) and how the scattering and spectral index affect what 
pulsars we are able to observe. Finally, in Section 5, we sum- 
marize our main conclusions. 



Number of beams 
Polarizations/beam 
Centre Frequency (MHz) 
Frequency channels 
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Galactic longitude range 
Galactic latitude range 
Number of survey pointing 
Sampling interval 
Observation time/pointing 
Limiting sensitivity 



^obs 



-60° to 30° 
-0.25° to 0.25° 
2560 
125 its 
1055 s 
0.15 mjy 




Angular Offset (degrees) 



Figure 1. Four pointings of the MMB receiver tessellated to give 
complete coverage over a whole patch of sky. One pointing is 
coloured black to show the pattern of the seven beams of the 
receiver, each of half-power width 0.053 degrees. 



2 SURVEY EQUIPMENT AND PROCESSING 
SYSTEMS 

2.1 Collecting Data 

This survey utilised the seven-beam Methanol Multibeam 
(MMB) receiver at a central frequency of 6.5 GHz mounted 
upon the 64-metre Parkes radio telescope. The seven beams 
of the MMB can be tessellated in order to cover the en- 
tire survey region, as shown in Fig. [I] with each beam 
having a half-power width of 0.053 degrees. The original 
plan was to survey a thin strip of the Galactic plane from 
-60° < I < 30°, |&| sC 0.5°. Subsequently, due to lack of 
success in finding previously unknown pulsars, the latitude 
range was reduced to |6| ^ 0.25°. Data were collected from 
February 2006 to September 2007. Table [l] lists the survey 
parameters and Fig. [2] shows the area of sky covered by the 
survey observations, with each seven-beam pointing marked 
by a single cross. For \b\ ^ 0.25° the survey is 85% complete 
and for 0.25° < |6| ^ 0.5° the completed fraction is 31%. In 
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Galactic Longitude (degrees) 



Figure 2. All survey pointings for the Galactic-plane portion 
of the MMB survey plotted in Galactic coordinates. The dense 
coverage in the range |6| ^ 0.25° contrasts with that for 0.25° < 
|6| ^ 0.5° where observations were not completed. 



addition, as shown in Fig. [3] a small area of ±0.3° around the 
Galactic centre was surveyed in July 2006 over 16 separate 
observations, with an integration time of 16,800 seconds per 
pointing, providing approximately four times the sensitivity 
of the observations of the Galactic plane region. Tables of 
the central beam positions in each survey are included in 
the supporting online material. 

A large observing bandwidth of 576 MHz, covered by 
a series of 192 filterbank channels each of width 3 MHz, 
was used in order to counteract the effects of interstellar 
dispersion. This dispersion causes a time delay in seconds, 

8.3 x 10 3 DM^" 3 A^ 



At : 



(3) 



across an observing bandwidth of Av which is centred at 
frequency v (both in MHz). During each pointing, the data 
were one-bit sampled every 125 fis, written to a local com- 
puter disk at the Parkes site and, once the pointing was 
completed, were written to DLT-S4 tape. This amounted to 
~5 TB of data, which were subsequently processed offline 
at the Jodrell Bank Observatory, while a second copy of the 
survey was kept on tape at the the Australia Telescope Na- 
tional Facility for backup and for quick access to data in the 
event of discoveries. 
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Figure 3. The position of all observations from the Galactic- 
centre portion of the MMB survey, from tessellating the MMB 
receiver in 16 different positions. 



in equation (JHjl, are properties of the pulsar to be detected, 
rather than survey parameters. 

Fig.[4]shows the computed sensitivity curves for disper- 
sion measures of and 1000 cm -3 pc respectively, using a 
mean sky temperature of 0.5 K, assuming S/Nmi n — 8. The 
intrinsic width of the pulse is assumed to be 5% and we use 
the scattering model of |Bhat et al. ( 2004 1 , shown in equa- 
tion |6j. The curves indicate that at periods greater than 
about 10 ms, the sensitivity is approximately 0.15 mjy. It is 
also evident that the large change in DM between the two 
curves results in a very small change in limiting sensitivity 
at periods above ~100 ms, which is an effect of the high 
observing frequency used in this survey, discussed further in 
Section 13.11 

For the Galactic-centre region, the sensitivity is affected 
by the background temperature of the Galactic centre itself. 
Using the maps of Seir adakis et al.| ( |1989| ), we estimate this 
contribution to be ~4 K in the outer regions of the survey, 
~12 K in the inner regions and ~90 K at the GC. The 
detection threshold of S/N^in — 9 is then ~28 /xJy (outer 
regions), ~32 nJy (inner regions) and ~81 /iJy at the GC. 



2.2 Survey Sensitivity 

The limiting flux density, in mjy, of a pulsar search (see 
|Dewey et aL]p85] | can be calculated by 

_ (g/JV) min /?T sya / W eff 
min ~ G^n p t ohs Au VP-W eff - W 

In this expression (S'/A r ) m i n is the minimum detectable 
signal-to-noise ratio, j3 is the degradation factor which takes 
into account effects such as signal digitisation, T sys is the 
system temperature (defined as the sum of the sky, spillover 
and receiver temperatures), G is the telescope gain, n p is 
the number of polarisations summed (n p = 2), t h s is the 
observation time and Av is the bandwidth in MHz. P, the 
pulse period, and W c ff, the effective pulse width as defined 



2.3 Processing the data - the first pass 

The processing pipeline used the SIGPROC-4.2 software 
package^] to perform dedispersion of the data up to a DM 
of 2000 cm -3 pc — chosen to ensure the survey was capa- 
ble of observing pulsars at the Galactic centre — using a 
DM step size of 7.18 cm -3 pc, derived from equation |3j), en- 
suring that the residual dispersive time delay between the 
highest and lowest frequency channels is no greater than 
one sample. These dedispersed time series were fast Fourier- 
transformed (FFT), and the spectrum was searched for sig- 
nificant peaks, as were spectra with 2, 4, 8 and 16 harmonics 
summed in order to ensure power was obtained from narrow 

1 http: / / sigproc.sourceforge.net / 
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Tobs=17.6 minutes 



Period (ms) 

Figure 4. Minimum detectable flux density (in mjy) of a source 
as a function of the period (in ms) for a detection at the 8cr level 
in the MMB survey when placed at a DM of and 1000 cm' 3 pc, 
calculated using equation Q. The sensitivities were calculated 
for integration times of 17.6 minutes and 280 minutes, for the 
Galactic-plane and Galactic-centre surveys respectively. 



pulses (Manchester et al. 1995). The Pulsar Hunter package^] 
was then used to search in period and DM space around the 
values output from the FFT search, in order to maximise 
the S/N ratio of each candidate. 

The pipeline was run on 'DCore', a beowulf cluster with 
72 processors, and candidates were inspected with JReaper 
( Keith et al.|2009| , a graphical tool for selecting candidates 
based upon user-definable parameters. This first pass yielded 



only one new discovery, J1410-6132 ( |0'Brien et al.||2008[ ) 



discussed in Section 4.1. 

One of the possible reasons why the first pass through 
the data found only one pulsar may have been the lack of 
dispersion delay, making it difficult to distinguish astrophys- 
ical signals from RFI at this high frequency (as discussed in 
Section 3.1 1. This resulted in an overwhelming number of 



spurious candidates, making the task of finding good candi- 
dates which should be reobserved very difficult. Therefore, 
it was decided that a second pass through the data would 
be useful, after introducing some techniques to remove, or 
reduce the effects of, these RFI signals. 



2.4 Reprocessing the data 

In reprocessing the data, much of the pipeline was very sim- 
ilar to that used originally. However, the processing was car- 
ried out on 'Hydra', a 108-node cluster, with each node con- 
sisting of dual quad-core processors with 4GB of memory. In 
order to take advantage of this processing power, a threaded 
tree algorithm written by one of us (Bailes) was used, which 
performed dedispersion many times faster than was possible 
with 'DCore' and the original dedispersion code. 

When preparing the data for processing, we decided not 
to use the time-domain clipping algorithm that was previ- 
ously employed by this and other surveys, since the discovery 
that this algorithm introduced periodic signals to the data. 
These signals were apparent even at high DM, polluting the 
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Figure 5. The normalised signal-to-noise ratio of a periodic sig- 
nal, with a duty cycle of 5% and period 500 ms, as the difference 
between the true DM value and the trial DM value varies. This 
is shown for both the MMB survey (dashed line) and the PMPS 
(solid line). The DM step sizes shown are those used for process- 
ing data in the respective surveys. 



search output. By removing this step, the effects of RFI at 
low DM were increased, however the sensitivity to these ar- 
tificial signals was removed. Steps were taken to reduce our 
sensitivity to RFI in other ways, outlined in the following 
sections. 

Instead of making changes to the processing pipeline, 
alterations were made in the post-processing stage, as out- 
lined in Section [3^2] For example, the data were not searched 
in acceleration space to look for pulsars in short-period bi- 
nary systems. The additional computation required by such 
a search is large, while this technique would also produce 
enormous numbers of candidates, something that we were 
attempting to avoid with this reprocessing. 



3 SEARCHING FOR PULSARS AT HIGH 
FREQUENCIES 

3.1 Dispersion as a discriminatory tool 

Interstellar dispersion, described by equation (J3j , is often 
a good way to distinguish between extra-terrestrial sources 
and terrestrial sources (that is, RFI), since one would ordi- 
narily expect that any periodic sources of RFI would not fol- 
low the dispersion law, and that they would peak in strength 
at a DM of cm -3 pc. At high frequencies, however, where 
the effects of dispersion are small, a signal with zero disper- 
sion will reduce in S/N very slowly as trial DM increases - 
the simulated pulse in Fig. [5] falls in S/N by less then one 
percent when the difference between the true DM and the 
trial DM is 50 cm -3 pc. By comparison, at an observing fre- 
quency of 1.4 GHz (in, for example, the PMPS), the same 
simulated pulse falls by ten percent in S/N when the true 
and trial DMs differ by only 4 cm -3 pc. Due to their short 
periods, this effect should not reduce sensitivity to millisec- 
ond pulsars (MSPs), for whom the dependance of S/N with 
DM is very strong. 

The effect of the small amount of dispersion delay at 
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Figure 6. The normalised signal-to-noise ratio of a signal with 
period 100 ms, for trial DM values every 25 cm -3 pc, at the MMB 
observing frequency of 6.5 GHz. The dotted, dashed and solid 
lines show S/N for signals with W/P 50%, 10% and 5% respec- 
tively. 



Fig.[6]shows S/N against trial DM error for signals with 
duty cycles of 50%, 10% and 5%; it can be seen that S/N 
drops off more quickly for narrower pulses. Thus, a filter 
was applied to the candidates to reject any with a width 
greater than 10% of the period, which eliminated sinusoidal 
signals often indicative of RFI, which are often of high S/N 
for a very large range of trial DMs. In total, these two fil- 
ters removed ~15% of the candidates from the processing 
pipeline. 

During the inspection of the candidates, a further cut 
was made, removing all candidates with S/N < 9. This 
removed a further 50% of the candidates, with many of 
those remaining clustered around common RFI frequencies. 
It is common, during the inspection process, to visualise 
the candidates in Period-DM space (for an example using 
the JReaper software, see Fig. 2 in Keith et al. ( 2009 1), so 
that any 'clustering' of the candidates at RFI frequencies is 
highly apparent, making it easier for such candidates to be 
ignored. This reduced the number of candidates to ~ 10000, 
and given that it takes only a few seconds to classify a can- 
didate, this allowed visual inspection of this subset of the 
survey output. 



high frequencies is that the variation of S/N with trial DM 
value is low compared to observations at lower frequencies. 
It is normal to use a plot of the variation of S/N with DM 
to help identify good candidates, but for this survey these 
plots were of little value, except for the highest DM sources 



(see Section 4.1.11. In lower DM cases, one has to judge 
instead based purely upon whether the periodic signal is 
persistent in time and frequency, and upon the shape of the 
pulse profile, with little or no indication of whether a source 
is actually terrestrial in origin. This led to there being an 
overwhelming number of candidates, especially at low DM. 
A typical single beam from the survey would generate ~170 
candidates, leading to a total of over 3.5 million candidates 
in the survey. 



3.2 Reducing the number of candidates 

This large number of candidates is impossible to filter by 
eye and so some automated filters were designed to reduce 
the number to be viewed. As the MMB receiver has multiple 
beams, it allows us to compare each of the beams from the 
same observation, and assume that if a periodic signal is 
detected in many beams, then it is likely to be ground-based 
RFI. Therefore, any periodicities which appeared in two or 
more beams from the same pointing and had a period which 
matched within a tolerance of 10 ns were rejected. 

As shown in equation Q, the S/N of a periodic signal 
decreases as the effective pulse width increases for a given 
set of observing parameters. The effective width is given by 



+ At 2 + t 2 



(5) 



where W sam p is equal to the sampling rate and r s is a distor- 
tion in the pulse width due to scattering in the ISM (given 
by equation |6|). Equation Q can be used to give At, where 
Av is now the bandwidth of a single filterbank channel. This 
broadening of the pulse arises due to the dedispersion pro- 
cess and the finite width of the filterbank channels in the 
observing system. 



4 RESULTS & DISCUSSION 
4.1 New Discoveries 

4-1.1 Galactic Plane Region 



The first processing of the data, by O'Brien et al. (2008), 



uncovered one previously unknown pulsar, the young and 
highly energetic PSR J1410-6132. |Q'Brien et~aL] p008| 
noted that the pulsar was within the error box of the 
EGRET source 3EG J1410-614 7 and it is within that of 
the Fermi so urce J1410.3-6128c (The Fermi-LAT Collabo^| 
|ration||2010|. The pulsar is regularly observed at Parkes as 



part of the timing support for the Fermi mission ( Wcltevrede 
et al.|[2010| ), but as yet pulsations in gamma-rays have not 
been detected ( |Abdo et al.|2010 |. 

One further pulsar, PSR J1834-0812, with a period of 
491 ms and a DM of 1023 cm~ 3 pc was discovered during the 
reprocessing of the Galactic-plane data (detailed in Section 
2.4 1. As seen in the top panel of Fig. [7] the pulse profile 
at 6.5 GHz is narrow with a duty cycle of only 2%, which 
ensured that the S/N was significantly lower away from the 
true value of DM, as shown in Fig. [8] 

The pulsar has been observed at a frequency of 4.85 GHz 
using the Effelsberg 100-m radio telescope, and at 1.4 GHz 
with the Lovell telescope. The profile is extremely scat- 
ter broadened at this frequency, with a scattering time of 
~150 ms. The pulse profiles at 4.85 and 1.4 GHz are shown 
in the bottom two panels of Fig. [7] 

4-1.2 Galactic Centre Region 

Processing of the data taken around the Galactic centre was 
completed in August 2006 and resulted in an excellent can- 
didate. The candidate was subsequently confirmed to be a 
pulsar, PSR J1746-2850, with a period of 1077 ms and a DM 
of 963 cm -3 pc in an observation made on October 29, 2006. 
Subsequently, this pulsar was also seen by |Deneva et al.| 
(2009), who have published a full timing solution. A folded 
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Table 2. Measured and derived parameters for each of the pulsars discovered in the survey. Numbers in parentheses signify the error in 
the last digit. Values for J1746— 2850 are taken from the timing solution of |Deneva et al.] ( |2009| ) , while those for J1834— 0812 are based 
on an interim timing solution spanning ~150 days. 



Pulsar Name 




J1746-2850 


J1834-0812 


Right Ascension (J2000) 




17:46:06.6(2) 


18:34:29.8(9) 


Declination (J2000) 




-28:50:42(5) 


-08:12:00(100) 


Period (s) 




1.0771014910(4) 


0.4911415563(2) 


Period Derivative 




1.34311(2) xl0~ 12 


9.96(7) xl0~ 15 


Dispersion Measure (cm" 


- 3 pc) 


962.7(7) 


1020(50) 


Flux Density at 6.5 GHz 


(mjy) 


0.14(3) 


0.10(2) 


Spectral Index 




-0.3 


-0.7 


DM Distance (kpc) 






9.6(1.2) 





1000 1200 1400 

Trial DM (cm" 3 pc) 



0.4 0.6 
Pulse Phase 



Figure 8. S/N of PSR J1834-0812 for different trial DM values 
used in the data processing pipeline. Note that the processing 
pipeline discovery slightly over-estimated the DM, hence the peak 
is positioned here at 1175 cm - 3 pc 



Figure 7. Pulse profiles of PSR J1834-0812 at observing fre- 
quencies of 6.5 GHz (top panel), 4.85 GHz (middle panel) and 
1.4 GHz (bottom panel), with the peak of each profile positioned 
at phase 0.1. The increase in pulse width as observing frequency 
is reduced, due to scattering, is clearly visible, particularly at 
1.4 GHz. 



pulse profile of this pulsar from a confirmation observation 
taken with the MMB receiver is shown in Fig. [9] 

Parameters for PSRs J1746-2850 and J1834-0812 are 
shown in Table [2] including DM distances estimated using 
the Galactic electron distribution model of Cordcs & Laziol 
(2002). As the parameters for PSR J1834-0812 are based 
on timing data that does not yet span one year, increased 
precision will be obtained with further observations and as 
period derivative and position become non-covariant. 




Pulse Phase 



Figure 9. Pulse profile of PSR J1746-2850, discovered in the 
Galactic-centre portion of the MMB pulsar survey. 



4.2 Known Pulsars 

Within the survey region, there are 113 previously known 
pulsars. Of these 113, 16 were detected. 



4-2.1 Non-detections 

A total of 62 pulsars fall below our detection threshold (tak- 
ing the offset from the beam centre into account) for any sen- 
sible value of the spectral index (i.e. where a < —0.5). For 
the remaining 16 pulsars, we can compute an upper bound 
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t he flux density at 1.4 GH z, as 
( Manchester et al.|2005| . An 



on the spectral index by using 
given in the pulsar catalogue! 

upper limit to the 6.5 GHz flux density was calculated as- 
suming a sensitivity of 0.2 mjy at the beam centre and the 
half power beam width of 0.053 degrees. These calculated 
values are shown in Table |3] 

Of the limits set on the spectral index, the steepest is set 
at -2.0, for PSR J1809-1943. This would not be an unusu- 
ally steep value for the spectral index, well within the esti- 



mated distributions of both Maron et al. ( 2000 \ and Lorimer 
|et al.| ( [T995| . Therefore, we can be confident that the flux 
densities are too low to have made them visible in the survey 
data. 



4.2.2 Detections 

Only one of the pulsars, PSR J1746— 2856, was undetected 
in the processing but later found when folded. This is likely 
the result of strong RFI during that observation. 

For those pulsars that were detected, Table [4] shows 
the S/N from the search pipeline (S/N Search), and that 
obtained when the data were folded at the period of the 
pulsar (S/N Fold). The other columns show the flux density 
at 1.4 GHz dLorimer et al. 



1995J, 4.85 GHz (obtained during 



the work of |Maron et al.| ( 2000[ )) and as estimated from our 
survey detections, using equation Q. 

We find that there is a large range of values for the 
spectral index, between 0.03 and —2.30 with an average of 
—0.94. This value is shallower than that found by |Maron| 



et al. (20001 who investigated pulsar spectra between 0.4 
and 4.9 GHz. However, by observing at a high frequency for 
a fixed duration, we are biased to observing pulsars which 
have a flat spectral index, or that are sufficiently bright to 
be observable at this high frequency. Since the majority of 
pulsars have flux densities of a few mjy or less at 1.4 GHz, 
it is not unexpected that our sample should contain many 
pulsars with a shallow spectral index. 



4.2.3 PSR J1550-5418 

PSR J1550— 5418 is a radio-emitting magnetar ( Camilo et al. 



20071. The pointing containing this source was taken on 
April 23rd 2007, 46 days prior to the observation by | Camilo] 
et^lll((2008l in which PSR J1550-5418 was discovered. Un- 



fortunately, the data were not passed through the pipeline 
until after the announcement of the discovery, in the radio 



at 1.4 GHz, by Camilo et al. (2007 k The flux density cal- 



culated from our observation does, however, agree with the 
observations of Camilo et al. (2008), given both the lack of 



a flux density calibrator for this single observation and the 
time variation of the flux density for that source. Since our 
observation fell between the initial observations, in the X- 



ray, of Gelfand & Gaensler ( 2007 1 and the radio detection of 



Camilo et al. ( 2008 1, we are only able to conclude that radio 
emission commenced at least 46 days before the published 
discovery. 



http: / / www.atnf.csiro.au/research / pulsar / psreat 



Table 3. Basic parameters and flux densities at 1.4 GHz for all 
undetected pulsars in the survey region with a limit on the spec- 
tral index steeper than —0.5. The last column, ai; m , gives an 
upper limit for the spectral index. It is calculated assuming a 
maximum flux density at 6.5 GHz of 0.2 mjy for a source at the 
centre of the survey beam. 



Pulsar Name 
(J2000) 



Period 
(seconds) 



DM 
(cm -3 pc) 



<Sl400 
(mjy) 



"lh 



J1512- 


-5759 





128694 


628.70 


6.0(6) 


-1.7(1) 


J1614- 


-5048 





231694 


582.80 


2.4(3) 


-1.5(1) 


J1648- 


-4458 





629632 


925.00 


0.55(7) 


-0.6(1) 


J1705- 


-4108 





861067 


1077.00 


1.3(2) 


-1.0(1) 


J1717- 


-3737 





682419 


525.80 


0.69(8) 


-0.76(13) 


J1722- 


-3632 





399183 


416.20 


1.60(17) 


-1.1(1) 


J1730- 


-3350 





139460 


259.00 


3.2(3) 


-1.8(1) 


J1739- 


3131 





529441 


600.10 


4.9(5) 


-1.4(1) 


J1740- 


-3052 





570310 


740.90 


0.7(2) 


-0.8(3) 


J1741- 


-3016 


1 


893749 


382.00 


2.3(5) 


-0.6(2) 


J1756- 


-2435 





670480 


367.10 


2.0(2) 


-1.5(1) 


J1809- 


1943 


5 


540249 


178.00 


4.5(5) 


-2.0(1) 


J1818- 


1519 





939690 


845.00 


2.1(4) 


-1.4(2) 


J1818- 


-1541 





551134 


690.00 


1.0(2) 


-0.9(2) 


J1819- 


1510 





226539 


421.70 


0.6(1) 


-0.6(2) 


J1822- 


-1400 





214771 


651.10 


0.80(9) 


-0.6(1) 



4.3 Spectral index distribution 

We can check for the completeness of our survey in the fol- 
lowing way. According to |Maron et al.| ( |2000| ), the spectral 
index distribution is Gaussian, with a mean of —1.8 and a 



(1995) 



standard deviation of 0.2, whereas in Lorimer et al 
the mean is —1.6 and the standard deviation 0.3. We take 
each of the 113 previously-known pulsars in our survey re- 



gion (gathered using the pulsar catalogue, Manchester et al. 
( |2005[ )) with a known flux densities at 1.4 GHz (leaving 103 
in the sample), and compute the probability of their detec- 
tion given the Maron and Lorimer distributions. 

To do this we choose a pulsar at random from the 103 
pulsars which remain, and assign it a spectral index cho- 
sen at random from a gaussian distribution with mean and 
standard deviation given by the Maron and Lorimer distri- 
butions. Using this spectral index, we then scale the flux 



density at 1.4 GHz (given by the pulsar catalogue, Manch- 



ester et al. ( 2005 1) to 6.5 GHz and infer whether this pulsar 
would be above the detection threshold of 0.2 mjy. Repeat- 
ing many times and normalising the results to 115 pulsars 
(since two have been discovered), we obtain an expected 
number of detected pulsars for each distribution. 

The expectation is that we should have found 12.6 
pulsars for the Maron et al. distribution and 16.1 for the 
Lorimer et al. distribution. The fact that we detected 18 
indicates that a) our survey is reasonably complete; and 
b) there is evidence for there being a small excess of pul- 
sars with flatter spectral indices than predicted by either 
Lorimer et al.| ( |1995[ ) or |Maron et "ah] |2000 ) , but no signifi- 
cant population with extremely shallow spectral indices. 
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Table 4. S/N values and flux density at 1.4, 4.85 and 6.5 GHz, and the computed spectral index (from 1.4 to 6.5 GHz) for the pulsars 
detected in the survey data. The final column shows the measured scattering timescale at 1.4 GHz and its error. 



Pulsar Name 


Period 




0/ IN 


^ /i\r 

Of IN 




S4.85 


06.5 


Q! 




(J2000) 


(seconds) 


(cm -3 pc) 


Fold 


Spectral 


(mJy) 


(mJy) 


(mJy) 




(ms) 


J lOl U U^Oa 




983.3 


41.5 


21 5 


0.74(8) 




U. 1 (L) 


— 0.0(1) 


29(8) 


TIQ07 fiOOO 




O-LO.O 


A 

oy .ft 


oy.u 












11 ^41 — (S22D 




717.3 


22.2 


13.5 






0.8(2) 


—0.5(1) 


8.3(3) 


J1406-6121 


0.213075 


542.3 


13.8 


11.1 


0.36(5) 


— 


0.20(5) 


-0.4(2) 


9(1) 


J1410-6132 


0.050052 


960.0 


12.4 


7.7 


6(1) 


— 


0.6(1) 


-1.5(2) 


33(15) 


11 531 — 561 n 


084202 


110.9 


15.6 


15.4 


0.60(7) 




v.'iyi-) 


—0.2(2) 




J1550— 5418 


2.069833 


830.0 


231 


199 


3.3(3) 




9(2) 


+0.7(1) 




J1644-4559 


0.455060 


478.8 


703 


175 


310 




12(2) 


-2.1(1) 


2.45(2) 


J1702-4128 


0.182136 


367.1 


25.6 


24.3 


1.10(12) 




0.6(1) 


-0.4(1) 


7.5(4) 


J1707-4053 


0.581017 


360.0 


10.7 


10.3 


7.2(7) 




0.19(5) 


-2.3(1) 


26.1(9) 


J1740-3015 


0.606784 


152.2 


18.5 


11.2 


6.4(7) 


1.08(0) 


1.7(3) 


-0.9(1) 




J1746-2856 


0.945224 


1168 


11.01 


N/A 






0.10(2) 






J1757-2421 


0.234101 


179.5 


15.1 


16.3 


3.9(4) 




1.0(2) 


-0.9(2) 




J1801-2304 


0.415796 


1074 


24.1 


16.8 


2.2(2) 


0.27(3) 


0.5(1) 


-1.0(1) 


99(9) 


J1803-2137 


0.133617 


234.0 


85.4 


81.7 


7.6(8) 


5.5(1.9) 


3.1(6) 


-0.6(1) 


5.5(7) 


J1809-1917 


0.082747 


197.1 


36.7 


35.7 


2.5(5) 




1.6(3) 


-0.3(2) 




J1832-0827 


0.647293 


300.9 


14.6 


13.0 


2.1(2) 


0.14(0) 


0.20(5) 


-1.7(1) 




J1834-0812 


0.491101 


1023 


18.0 


12.0 


0.3(1) 




0.11(2) 


-0.7(3) 


150(100) 



4.4 Scattering 

By taking the 115 known pulsars (including our discoveries) 
in the survey region, it is possible to analyse the impact 
of scattering by the ISM. This is especially important at 
low Galactic latitudes, where scattering is most severe and 



which was not well covered by the Bhat et al. ( 2004 ) survey. 



Scattering by the ISM causes a pulse profile to be convolved 
with an exponential tail of decay timescale r s , which was 
found to follow the empirical relation (for r s in milliseconds), 

log r s = -6.46 + 0.154 log DM + 1.07(log DM) 2 - 3.86 log v, 

(6) 

at a frequency of v GHz, for a sour ce with dispersion mea- 
sure DMcm" 3 pc ( |Bhat et al.|2004 >. 

In order to measure t s we used the pulse profiles at 
6.5 GHz obtained from the survey and, where available, 
profiles at 1.4 GHz obtained elsewhere. This amounted to 
12 pulsars in total. We assume that the scattering in the 
6.5 GHz profile is negligible compared to that at 1.4 GHz 
and convolve this 6.5 GHz profile with an exponential scat- 
tering tail of decay time t 3 , fitting against the 1.4 GHz pro- 
file. We ignored values of r s less than 2 ms where the results 
are not reliable, and further ignored any pulsars in which 
there was noticeable pulse profile evolution between the two 
frequencies. The values of r 3 , and their standard deviation, 
for the remaining 10 are given in column 10 of Table [3] 
Fig. [To] shows DM vs r s for this group of pulsars, with the 
relationship given by equation |6]| plotted. 

Fig. [TO] shows that the scattering tail size is clustered 
in the region from 10-100 ms, with the lower boundary im- 
posed, as previously stated, by insufficient time resolution 
in the data to measure any small difference in profile shape 
at lower DMs. The higher boundary comes from the lack of 
known sources at DM > 1000 cm -3 pc, something this sur- 
vey hoped to address when it began. What it does show, 
however, is that high-DM (~ 1000 cm -3 pc) sources would 



700 800 900 1000 



DM (cm^-S pc) 



Figure 10. Characteristic width of an exponential scattering tail 
convolved with a pulse profile at 6.5 GHz in order to produce the 
best-fit to the pulse profile at 1.4 GHz. The dotted line shows the 
model of |Bhat et al.| |2004) for a frequency of 1.4 GHz. 



be so highly scattered at 1.4 GHz that it is unlikely the 
PMPS would have discovered them, no matter how high 
their flux density. The fact that so few were found in this 
survey, indicates that the amount of obscuring material be- 
tween ourselves and these putative sources, coupled with 
the lower flux density, prevented their discovery. The four 
highest-DM sources in Fig. [To] also he below the model of 



Bhat et al. (20041. Clearly, when observing at low frequen- 



cies, surveys are biased to detecting pulsars with less scat- 



tering, and thus the model of Bhat et al. ( 2004 1 can be seen 



as a lower limit to scattering. However, our results suggest 
that the relation between scattering and DM may be flatter 
than previously thought. 
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Table 5. Parameters used to simulate the pulsar population used 
in Section l4.4.1l 

Galactic latitude scale height 330 pc 

Minimum luminosity O.lmjykpc 2 

Maximum luminosity 10000 mjykpc 2 

d log iV/d log L -0.8 

log (Mean of P distribution (ms) ) 2.4 

log (Std. dev. of P distribution (ms) ) 0.34 

Number of pulsars detectable by PMPS 1033 



4-4-1 Simulated Population 

To fully understand the interplay between the spectral in- 
dices and degree of scattering, we produced a simulated pop- 
ulation of pulsars in which both of these factors could be 
allowed to vary. Using the PSRPOP software^ for each re- 
alisation, a simulated population of pulsars was produced 
with a different mean and standard deviation of the spec- 
tral index. All other variables were held fixed at the values 
given by |Lorimer et al.| ( |2006[ ) , and are listed in Table [5] 
The population was increased in size until it contained a 
subset of 1033 pulsars with flux density and position in the 
Galactic plane such that they would be detectable by the 
PMPS. As the PMPS is one of the most-studied pulsar sur- 
veys, responsible for the discovery of a large fraction of the 
pulsar population, it is suitable for constraining the popula- 
tion of pulsars in the Galactic plane. The survey parameters 
are outlined in Table [6] PSRPOP models the population of 
'normal' pulsars, and so our simulation is unaffected by the 
acceleration effects which are most commonly observed for 
MSPs. The acceleration of a pulsar in its orbit can smear the 
pulse out such that it is not discovered in a blind search, if 
the data is not searched in acceleration space. 

PSRPOP produces values of the flux density at 1.4 GHz, 
the pulse width and distance for each pulsar, based upon 
the analysis of the known pulsar population presented in 



Lorimer et al. ( 2006 1. Using these, the flux density and pulse 



width were then recalculated for the entire simulated popu- 
lation at 6.5 GHz, each time assuming a different power-law 
dependance of the scattering tail size with frequency, from 
equation |6| - labelled as "scattering index" in Figs. |ll| and 
|12| From these values, the number of pulsars which would 
be visible in the MMB survey region could be calculated for 
each simulated population. 

We could then find all simulated populations in which 
the number of detected pulsars is within one standard devia- 
tion of the number that we detected. These results are plot- 
ted as histograms in Fig. |ll| and on a greyscale in Fig. |12| for 
each pair of simulation parameters. From Fig.[l2]we can ex- 
clude certain regions of the parameter space that are clearly 
inconsistent with our results. 

For example, there is a clear favouring of spectral in- 
dices of ~ —1.5 in the histogram, as we might have expected 
given the results found in Section [4.3| It is, however, diffi- 
cult to constrain the scattering index in any way, as some 
simulations fit our observations, regardless of the value it 
takes. 

We can also put quantitative limits on the number of 
4 http://psrpop.sourceforge.net 



Table 6. Observational parameters for the Parkes Multibeam 
Pulsar Survey (PMPS) 

Number of beams 13 

Polarizations/beam 2 

Centre Frequency (MHz) 1374 

Frequency channels 96 X 3 MHz 

System temperature (K) 21 

Galactic longitude range 260° to 50° 

Galactic latitude range \b\ ^ 5° 

Sampling interval 250 fis 

Observation time/pointing 2100 s 

Limiting Sensitivity for centre beam 0.14 mjy 



pulsars with flat spectral indices, based on the small num- 
ber of detections in our survey, by simulating a population 
with a fixed spectral index at —0.5, and 0.5. By comparing 
the number of detections such a population would give com- 
pared to what we have observed, we can estimate a limit on 
how many such pulsars could be in the Galactic pulsar popu- 
lation. The simulations show that the population of pulsars 
with spectral index greater than —0.5 cannot exceed 10% 
and cannot exceed 5% for a spectral index of 0. 

Our findings have relevance to future high-frequency 
surveys planned with the phase- 1 Square Kilometer Array 
(SKA). Using the specifications as described in Schilizzi 



et al. (20071, we repeated this simulation for the case of 620 



15-metre dishes with system temperature of 35 K and effi- 
ciency of 65%, at an observing frequency of 10 GHz with a 
bandwidth of 500 MHz. In order to account for the apparent 
deviation of our results from those of iBhat et al. (20041, we 



obtained results using the scattering law described in equa- 
tion Ipl, and also using a linear fit to the data in Fig. 10 

Our simulation suggests that a survey of the inner 
Galaxy (|/| 45°, \b\ ^ 5°), using an integration time per 
survey observation of 1800 s, could expect to find ~ 2700 
pulsars, assuming a spectral index distribution with mean 
— 1.7, standard deviation 0.3 and a scattering index of —3.8. 
The predicted number of detections, n, changes by less than 
^Jn whether using our linear fit, or that of Bhat et al. ( 2004 1 , 
to model the dependance of scattering timescale with DM. 
However, we are yet to discover the true Galactic-centre pul- 
sar population, and it is likely that when we do so, the degree 
by which this population is scattered will be much larger, 
requiring the use of high observing frequencies. 

In order to compare the merit of such a survey to one at 
a lower frequency, a similar simulation was also performed 
for an observing frequency of 1.4 GHz. This survey would 
detect many more pulsars over same region, and would have 
a vastly improved survey speed. However, for Galactic lat- 
itudes \b\ ^ 2°, the high-frequency survey is predicted to 
find three times as many pulsars with DM > 900 cm -3 pc as 
the low-frequency survey. Hence these two survey types are 
complementary to one another, provided the high-frequency 
survey is targeted at discovering the Galactic centre popu- 
lation of pulsars. This is not possible with a low- frequency 
survey of the same region, due to the large amount of scat- 
tering introduced by the material between ourselves and the 
Galactic centre. 



10 S. D. Bates et al. 




Scattering Index Spectral Index Spectral Index Width 



Figure 11. Histograms showing the relative number of simulations (linear scale) consistent with our observations, for each of the three 
parameters. 




Spectral Index, a Spectral Index Width Spectral Index Width 



Figure 12. Grey scale plots showing the number of simulations that were found to be consistent with our results, for different combinations 
of simulation parameters. 



5 CONCLUSIONS 

A survey of the Galactic plane in the region —60° ^ / ^ 
30°, |6| ^ 0.25° was carried out between February 2006 
and September 2007 using the 64-metre Parkes radio tele- 
scope at a frequency of 6.5 GHz. The seven beam Methanol 
Multibeam receiver was used to increase the amount of 
sky coverage in a single pointing, and reduce the length of 
time required to complete the survey. Three pulsars have 
been discovered from two passes through the survey data; 
PSR J1410— 6132, a young pulsar with an extremely high 
DM of 960cm" 3 pc, PSR J1834-0812, also with a high 
DM of 1023 cm" 3 pc, and PSR J1746-2850, close in an- 
gular distance to the Galactic centre and with a DM of 
~ 1000 cm -3 pc. This small number of detections indicates 
that any population of pulsars with fiat spectral indices in 
the inner-Galaxy must be small. 

Analysis of the known pulsars in the survey region has 
shown that the decrease in intrinsic flux density of pulsars 
as a result of the spectral index removes any benefit ob- 
tained from the reduced scattering at this frequency. Of the 
known pulsars seen in the survey data, the majority have a 
shallow spectral index, and the few that do not are intrinsi- 
cally bright sources, which leaves them above our detection 
threshold, while one, J1550— 5418, a magnetar, has emis- 
sion with a negative spectral index, making it brighter at 
6.5 GHz than at the typical survey frequencies of 1.4 GHz 



or lower. The date of this observation places an important 
limit on when the source first emitted radio pulsations, since 



Gelfand & Gaensler ( 2007 1 were unable to observe pulsations 



in August 2006. 



The discovery of three pulsars with DMs in excess of 
900 cm -3 pc indicates that observing at higher frequencies 
does, indeed, allow one to probe at great distances in the di- 
rection of the Galactic centre. However, the reduction in dis- 
persive delay makes distinguishing which sources are likely 
to be astrophysical extremely challenging. We do find, how- 
ever, that there is some evidence that the relationship be- 
tween scattering timescale and DM is flatter than suggested 



by the model of Bhat et al. ( 2004 1. This strengthens the case 
for future surveys at lower frequencies, for example a new 
1.4 GHz multibeam survey at Parkes ( Keith et al.|2010 l and 
the all-sky survey using LOFAR (van Leeuwen & Stappers 



2009) 



Future high-frequency surveys of the Galactic centre re- 
gion with the SKA ( Smit s et al.|2009| could potentially dis- 
cover thousands of Galactic-centre pulsars. However, this 
will require a large observing bandwidth to increase the dis- 
persive delay across the band, and location of the telescope 
in a region with improved RFI conditions. 
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